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ABSTRACT 

We have identified seven (field) elliptical galaxies acting as strong gravitational 
lenses and have used them to measure cosmological parameters. To find the most likely 
value for (= 0,matter) and A, we have used the combined probabilities of these 
lens systems having the observed critical radii (or image deflection) for the measured 
or estimated values of lens redshifts, source redshifts, and lens magnitudes. Our 
measurement gives A = 0.641q;26 if + A = 1, and the flm = 1 model is excluded 
at the 97 % confidence level. We also find, at the 68 % (0 = 0)- 82 % {Q = 0.3) 
confidence level, that an open universe is less likely than a flat universe with non-zero 
A. Except for the possibility of strong perturbations due to cluster potentials and the 
systematic overestimate of the lens magnitudes, other possible systematic errors do not 
seem to influence our results strongly: correction of possible systematic errors seems 
to increase the significance of the result in favor of a non-zero A model. 

Subject headings: cosmology: theory - cosmology: observations - cosmology: 
gravitational lensing - galaxies: evolution 
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1. Introduction 

Recently, non-zero cosmological constant (A) models have found increased popularity (e.g., 
Ostriker & Steinhardt 1995) owing to the problem of the age discrepancy implied by the latest 
Hubble Space Telescope (HST) measurements of the Hubble constant: the ages of globular clusters 
are apparently larger than the age of the universe predicted by the standard Q.m = 1 model which 
is favored by standard inflationary theory (e.g., Freedman et al. 1994; for a brief overview of these 
arguments, see Rees 1996). 

In order to measure A, it has been suggested that strong gravitational lenses might be used, 
i.e. isolated galaxies or clusters of galaxies for which the gravitational potential results in multiple 
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imaging of a background object (Pacz3mski & Gorski 1981; Alock k, Anderson 1986; Gott, Park, 
& Lee 1989). Following these suggestions, the use of the lens number counts (or the optical 
depth) was advocated since this is very sensitive to A (Turner 1990; Fukugita, Futamase, &: Kasai 
1990; Fukugita et al. 1992). Maoz & Rix (1993) and Kochanek (1995) have applied this method, 
obtaining upper limits of A ^ 0.7. 

Kochanek (1992) has also suggested the lens redshift method, which, compared with the 
method based on lens counts, requires less presumptions about the properties of lenses and sources, 
properties which might bias the lens counts considerably (Helbig &; Kayser 1995; Kochanek 1992; 
Fukugita & Peebles 1995). Taking into account the selection effects which were neglected in his 

early study (Kochanek 1992: for a discussion on this selection effect, see Helbig &: Kayser 1995), 
Kochanek (1995) finds A < 0.9 at 2 a with a peak at A = 0.4. However, the estimated value of A 
in Kochanek (1995) is sensitive to the detection threshold which is not well understood, and thus 
it can not be considered very seriously at the present stage. 

It has been recognized that the mean splitting of the lensed images alone is useful for studies 
of the dynamical properties of lens galaxies, but not for the measurement of A (Turner, Ostriker, 
&; Gott, hereafter TOG84; Fukugita et al. 1992). However, when the mean separation is used 
together with other information such as the lens redshift, the lens magnitude, and the velocity 
dispersion of the lens galaxy, then the mean separation does become sensitive to A (Paczynski &; 
Gorski 1981; Gott et al. 1989; Kochanek 1992; Miralda-Escude 1991). In this Letter, wo will try to 
measure A using a method which we call the "lens parameter method" , which is basically similar 
to those discussed in the above references. 



2. Lens Parameter Method 

The commonly observed parameters for gravitational lenses are the lens redshift (zl), the 
source redshift (zs), the mean deflection of the lensed object (or similarly the critical radius 
{6crit))i the lens magnitude {mi)-, and the source magnitude {ms)- For some systems, one or two 
of these observational parameters may be missing. 

How sensitive is Ocrit to and A for a given set of zl, zs, and m^? To calculate Ocrit, we 
will adopt the singular isothermal sphere (SIS) model for the lens, along with the filled-beam 

approximation (see section 4), and the Faber- Jackson relation (Faber & Jackson 1976). The 
Faber- Jackson relation relates the velocity dispersion (a) of E/SO galaxies to their luminosities 
(L): a = cj*(L/L*)^. Here, we adopt a* = 225km/sec, j3 = 0.25 and M^Bt = —19.9 + 5log{h), 
taken from Kochanek (1992,1995). 

Then 9crit{zL, zsjItil) can be expressed as; 

Ocrit = 47r(c7./c)2 d{zs, ZL){d{0, Zl){1 + ZLffP/d{Q, Zg) lQ-^mm-M.-K{z)-Eiz)-2^) (-^^ 
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where d{zi,Z2) is the angular size distance between the redshifts Zi and Z2 in Mpc, m is the total 
apparent magnitude of the lens galaxy, K(z) and E(z) are the K-correction and the evolutionary 
correction, respectively, for the lens galaxy. The (E+K) correction is important, and to calculate it 
we will use the 1 Gyr burst model of Bruzual &: Chariot (1993) at the formation redshift Zfor = 10. 
This (E+K) correction is consistent with the results from the HST Medium Deep Survey (MDS) 
on the evolution of the luminosity function of elliptical galaxies (Im et al. 1996), which shows a 
brightening in luminosity by about 1 magnitude looking back to z ~ 1. 
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Figure 1: The 6crit-zs relation for two gravitational lenses HST14176+5226 and HST 12531-2914, assuming two 
different sots of cosmological parameters, i.e., (fim = 1,A = 0) and (fi^ = 0.2, A = 0.8). Observed values for Ocrit 
are indicated with the dsished lines. 

Fig. 1 shows the Oc-it ^ zs relation for the strong gravitational lens systems HST 12531-2914 
and HST14176+5226, taking parameters from Ratnatunga et al. (1995). The two curves show the 
model predictions under the adoption of different cosmological parameters, and the horizontal line 
shows the observed value (the source redshift is unknown). The value of Ocru is quite sensitive to 
A when values are known for zl, zg, and mi,. However, the uncertainty in the prediction is about 
a factor of 10°'^^, which arises mainly from the uncertainties in the Faber- Jackson relation and in 
the apparent lens magnitude. Hence, a single lens system such as HST 12531-2914 cannot be used 
alone to measure A. In order to set a useful limit on A with this method, a sample of at least five 
lenses is required (e.g., sec Kochanek 1992). 

In order to combine the information on cosmological parameters from all available lenses, we 
therefore construct a likelihood function which is the product of the probability of each lens having 
the observed value of OcHt for the given values of zg, zi, niL and the cosmological parameters. 
This probability Pi{9o; Zl, zs, mi, ri^j ^) is defined as, 

Pi{Oo; ZL, zs,mL,^m, ^) ~ / G{9o,9crit{z,zs,m{z)),a0)G{z,ZL,(JzL)G{m{z),mL,(Tjni^)dz (2) 
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where is the dispersion in the predicted logio{9crit) due to the uncertainty arising from the 
Faber-Jackson relation, together with other minor uncertainties, cj^^ is the uncertainly in zl, and 
is the uncertainty in m^. G{x,XQ,dx) is the Gaussian function with the mean of xq and the 
dispersion of dx. We adopt ag = 0.14 (or in terms of magnitude, a ~ 0.7) which is a combination 
of the uncertainties arising from the Faber Jackson relation (0.65 x 2 = 0.13: de Zeeuw &; Franx 
1991), M* (0.3 X 0.2 = 0.06: Marzke et al. 1994; Loveday et al. 1992), and the E+K correction 
(0.5 X 0.2 = 0.1: Im et al. 1996). When zs is not available (HST12531-2914), we also integrate 
Eq.(2) over zs, assuming a uniform distribution in redshift space. 

Finally, the likelihood function can be written 

= n Pj,normiOj ; ZL,j , Zs,j , mij ) (3) 
j 

where Pj^norm is the normalized probability of Eq.(2). 

We did not adopt the (3/2)°-^ factor (hereafter TOG factor) which was suggested by TOG84 
in order to account for the possible difference between the velocity dispersion of the underlying 
dark matter and the luminous material. Recent studies show that this factor is not necessary 
(Kochanek 1993,1994; Breimer Sz Sanders 1993; Franx 1993). Independently, we also checked the 
necessity of the TOG factor by considering the mean image splittings (see Section 4). 

The advantage of this method over the previous methods is the explicit use of the lens 
magnitude and the E+K correction, of which the latter has been observationally constrained only 
recently (Im et al. 1996; Pahre, Djorgovski, &: de Carvalho 1996; Bender, Ziegler, & Bruzual 1996; 
Barrientos, Schade, & Lopez-Cruz 1996). These measurements enable us to make a reasonably 
good estimate of the dynamical properties of each lens galaxy. The probability of each individual 
lens having its unique configuration can then be calculated based on these individual properties, 
so that we do not have to use statistical measurements (e.g., the luminosity function) which 
may decrease the dependence on the cosmological parameters when they are averaged over large 
numbers of objects. Although our method is, in principle, not as sensitive to the value of A as is 
the lens number count method (e.g., Fukugita et al. 1992), the latter method is possibly subject 
to greater uncertainties (see section 4). Our method is slightly more susceptible to a small change 
in one of the input parameters, but, in common with the lens redshift method, we have a smaller 
number of parameters than the lens count method. In this respect, our method has an edge over 
the latter. In particular, the properties of lens galaxies at high redshift (z ^ 1.5) are highly 
uncertain. They could be dusty enough that the result from the lens counts might be biased 
against the non-zero A model (Fukugita & Peebles 1995). In contrast, the lens parameter method 
uses lensing galaxies which lie at z ^1 (see section 3), and the method is thus less affected by the 
unknown properties of high redshift galaxies. 
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3. Sample Selection 

Gravitational lenses are selected using the following criteria. 

1) the strong lensing must be caused by a single galaxy lens. For example, we do not include 
2016+112 in our sample since there are two lensing galaxies in this system. Also, we have excluded 
lens systems which are clearly influenced by strong perturbations due to cluster potentials (e.g., 
0957+561, B1422+231). 

2) it must be known that the lens galaxy is likely to be elliptical. For example, we do not 
include B0218+357 in our study since there is good evidence that the lensing galaxy is a spiral or 
a late- type galaxy (Patnaik et al. 1993). 

3) the apparent magnitude and the rcdshift of the lens galaxy must be known or estimated to 
reasonable accuracy. Accurate values for and zl are important for estimates of the dynamical 
properties of the lens galaxy. 

4) For lens candidates that do not have a measured value for zs, we select only those which 
show distinctive features such as rings or crosses. 

Wc find that there are seven strong gravitational lenses that meet these selection criteria in 
the published literature, including objects found in our HST surveys (Table 1). B1422+231 is 
excluded from this list because of the possible cluster perturbation as well as the ambiguity in the 
lens redshift [zl = 0.64 from Hammer et al. 1995 vs. zl — 0.4 from Impey et al. 1996). For 
MG0414+0534, there have been speculations that the source redshift of the system is z ~ 1.00 
(Burke 1990; Kochanek 1992), but it was later established to be zs = 2.64 (Lawrence et al. 
1995), suggesting that the z ~ 1 measurement pertains to the lens galaxy (Surdej & Soucail 
1994). We have analyzed the archived HST observations of this system and find a preliminary 
result of Rf675W — -^F8i4VK = 1-5 ± 0.2 for the lens galaxy (Ratnatunga et al. 1996), suggesting 
that Zl = 1.2 ± 0.4, consistent with the previous estimates of zl = I, and hence we will adopt 
Zl = 1.2 ± 0.4 for this system. Finally, we have subtracted a few tenths of a magnitude from 
some of the quoted lens magnitudes in the literature, in order to correct for the total apparent 
magnitude. When the uncertainty in the lens magnitude is not quoted in the relevant reference, 
errors of about 0.3 magnitude are assigned to these lens galaxies. 



4. Results and Discussion 

In Fig. 2, we present the relative likelihood of our measurement against for two cases 
which are of cosmological interest: i) Q„i + A = 1 and ii) A = 0. Both likelihood functions are 
normalized with the maximum likelihood of case (i), and direct comparison of cases (i) and (ii) is 
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possible using Fig.2. When a fiat universe is assumed (case (i)), we find that A = 0.641^;^^, and 

we exclude the = 1 niodel with 97 % confidence. Also, a universe with A ^ 0.9 is excluded 
at the 95 % confidence level. If A = is assumed (case (ii)), then $7^ ~ is favored. The 
difference in the likelihood function between a fiat universe with A = 0.64 and an open universe 
with < < 0.3 is about 0.5-1. Hence, a flat universe with non-zero A is favored over an open 
universe at 68 % - 82 % confidence. 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 2: The likelihood of our measurement against Qm for two model universes: i) Qm + A = 1 (solid line), and ii) 
A = (dotted line). The relative likelihoods are normalized with the maximum of the likelihood of (i). Confidence 
limits are indicated on the right-hand side of the plot, assuming the distribution for the likelihood. 

Our result is only marginally consistent with the previous estimate of A < 0.7 based on the 
lens counts which strongly favored the zero A flat universe (Kochanek 1995; Maoz Sz Rix 1993). 
To see what might have caused the disagreement between our result and the previous results, we 
have investigated the possible systematic errors in our analysis and these are listed below: 

1) the flUed-beam approximation vs. the empty beam approximation 

To relate redshift to distance, the filled beam approximation assumes that light rays propagate 
through smoothly averaged spacetime. In reality, spacetime is inhomogeneous, and therefore 
the filled beam approximation may not be correct (e.g., Fukugita et al. 1992). To see how our 
result could be affected by the filled beam approximation, the analysis was repeated adopting 
another extreme assumption, namely the empty beam approximation. We find that the latter 
approximation does not change our result significantly, but strengthens our finding slightly in 
favor of the non-zero A model. 

2) Singular isothermal sphere vs. softened isothermal sphere. 

We can also assume different mass models for the lens, rather than the SIS model. Recent 
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studies show that the SIS model may be too simple to adequately describe the mass of E/SOs 
(Laucr 1988; Krauss & White 1992), although the size of the core radius may be small enough 
to be negligible (Wallington & Narayan 1993). If the softened isothermal sphere is used, the 
predicted Ocru will be a bit smaller than the predicted Ocrit with the SIS model. In Fig.l, this 
means that the predicted lines need to be shifted down along the Ocrit axis, making the Q,m = 1 
fiat model more inconsistent with the prediction. The adoption of the softened isothermal sphere 
model will thus strengthen our result. 

3) Morphological misclassification 

In our analysis, wc have assumed that each lens galaxy is an E/SO. This assumption may be 
wrong, and to estimate the bias introduced by treating a spiral galaxy lens as an elliptical galaxy 
lens, we repeated our analysis with the inclusion of one known spiral lens system (B0218+357) 
treated as an E/SO lens. This caused the result to be strongly biased in favor of the A = model, 
because of the small predicted 9crit of the spiral lens system, a result similar to the issue discussed 
in (2). Thus, if one of the seven lenses we used was actually a spiral galaxy rather than an 
elliptical, then the correction of it would only strengthen our result. 

4) Wrong lens magnitude 

Because the lens galaxy is much fainter than the lensed object in some cases, there is a 
possibility that the lens magnitudes are not well determined. Systematic overestimate of the 
lens magnitudes by more than 0.6 magnitude can bias the result against the Qrn = 1 model. 
Recent HST observations provide clues as to the accuracy of the ground-based estimates of lens 
magnitudes. The preliminary result by Falco (1995) from the HST observation of 0142 — 100 
gives an aperture magnitude of Rf675W = 19.66 ± 0.01 for the lens galaxy, while the original 
measurement by Surdcj et al. (1987) is R=19. For MG0414+0534, our preliminary analysis of the 
HST observation shows / ~ 21.4 it 0.15 for the lens galaxy (Ratnatunga et al. 1996), agreeing 
with the previous estimate of / ~ 21.08 — 21.36 from the ground (Schechter & Moore 1993). On 
the other hand, Impey et al. (1996) find rriL = 21.5 ± 0.3 in V for the lens galaxy of B1422+231. 
The ground based estimate is r = 21.8 =b 0.6 for this object (Yee &; Ellingson 1994; Yee 1995). At 
z ~ 0.4, y — r ~ 1 for E/SO galaxies, thus the observed ground-based lens magnitude for this 
system disagrees with the HST observation by about 1 magnitude. These three examples may 
indicate that the early ground-based measurements are not very accurate. Errors seem to go both 
ways, and hence there may be no systematic overestimate of the lens magnitudes. But there are 
only seven lenses in our sample, and it may be premature to say that there are no systematic 
errors in the lens magnitudes. 

5) The TOG factor 

With the TOG factor included, wc find that the peak of the likelihood function shifts to 
the region where both $7^ and A are very small, where our test becomes quite insensitive to the 
cosmological parameters. However, many studies have shown in different ways that the TOG 
factor is not needed (Kochanek 1995 and refs therein). In order to confirm earlier findings, we 
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analyzed the mean image deflections of the known lens systems with the known source redshifts. 

Using criteria (1) and (2) described in section 3, we find that there are 11 lens systems available for 
this analysis (see table 1 in Keeton &: Kochanek 1995). For these systems, we calculated the ratio 
dcrit,ohs/ < Gcriti^s) >■ Since < Ocriti^s) > is fairly independent of the cosmological parameters 
(TOG84; Fukugita et al. 1992), the average of 11 6crit,obs/ < Ocrit{zs) > values will be about 1 
if the TOG factor is not necessary and about 1.5 if the TOG factor is appropriate. We find an 
average value of 1.0 it 0.1, confirming that the TOG factor is not necessary. In order to test the 
TOG factor independently, the study of strong lenses at low rcdshift {z < 0.1) might be fruitful, 
since their lens parameters are then insensitive to the cosmological parameters. An optical survey 
which covers a large fraction of sky (e.g., SDSS) should be able to find a statistically significant 
number (~ 100) of such lenses. 

6) Cluster Perturbation 

Since elliptical galaxies preferentially live in a cluster environment, the gravitational potential 
of the lens may include a cluster component. The strong cluster perturbation generally increases 
the mean image deflection, and hence we tried to exclude such lenses from our study (See section 
3-(l)). Nevertheless, we can not completely exclude the possibility that some of the lenses in our 
sample include a considerable amount of cluster perturbation. If that has happened, our result 
could be biased against the rt,ri = 1 universe. To understand the possible contribution to the 
image splitting from the cluster potential, detailed modeling of the lens systems is desired using 
high resolution images from the HST, or else radio observations. 

7) Source redshift for HST14176+5226 

Crampton et al. (1996) have recently published a tentative source redshift for the lens system 
IIST14176+5226. A strong emission line is found at 5324 A, along with a possible weak emission 
feature at 6822 A. The strong emission line is very likely to be Ly a at z = 3.39 if the weak 
emission feature at 6822 A is real, and the latter can then be identified as CIV 1549. If the 6822 
A feature is not real, then the source object could be located at a redshift lower than zs = 3.4. If 
zs < 3.4 for the HST14176+5226, then the predicted 6crit will be reduced. This would bring the 
peak of the likelihood function toward the large A value, strengthening the result in favor of the 
non-zero A model (see Fig. 1) . 

8) E+K correction 

The adopted E+K correction assumes a formation redshift of Zfoj- = 10 with a 1 Gyr burst 
of star formation. We find that the E+K correction is most sensitive to the value of Zfor, and 
insensitive to the other parameters. If we adopt Zfor > 10, then the result changes insignificantly 
towards the zero A model. When Zfor < 10, the result changes in favor of the non-zero A model, 
and the change is significant when Zfor < 2. If Zfor = 1-5, A could be as large as A ~ 0.8. 

If our result is an overestimate in the value of A, then there must have been a large systematic 
overestimate of the lens magnitudes and/or there are strong cluster perturbations. On the other 
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hand, if the lens counts have led to an underestimate in the value of A, then that could have 
been caused by: i) the dusty nature of high redshift elliptical galaxies (see section 2 for more 
discussion), ii) a decrease in the number density of ellipticals as a function of look-back time, as 
expected if most elliptical galaxies were created via major merging events (Im et al. 1996,1997; 
Baugh, Cole, k. Prenk 1996; Kauffmann, Chariot, k, White 1996), and/or iii) other uncertainties 
in the properties of lens galaxies, such as the LF and the dynamical properties of the low mass 
ellipticals. Future HST observations of faint galaxies, as well as the accumulating redshift data 
from ground based telescopes, will hopefully put stringent constraints on elliptical galaxy evolution 
at z > 1. These data will possibly give us indications as to why the results from the lens counts 
have strongly favored the zero A model while our result strongly rejects the flat universe with 
A = 0. It is noteworthy that neither method strongly rejects the low universe. 



5. Conclusions 

We have described and applied the lens parameter method to measure cosmological parameters 
using strong gravitational lenses. Using seven strong lenses each with an identified lens galaxy, we 
find that a model universe with A ~ 0.65 and low J7 is favored and that the flat model with A = 
is excluded at greater than 95 % confldence. A universe with low O and A = can be marginally 
excluded with respect to the flat universe with a non-zero A at 68 % — 82 % confldence. Our 
result is not biased in favor of a non-zero A model due to any conceivable systematic errors, except 
for possible strong perturbations from cluster potentials, and systematic overestimate of the lens 
magnitudes. Future HST observations should uncover new lens systems with measurable lens 
properties suitable for this kind of study, and they should also provide a better understanding of 
the known lens systems. We should therefore be able to get a stronger constraint on A in the near 
future. 

The HST Medium Deep Survey is funded by STScI grants G02684 et seqq.. We would like to 
thank the other members of the Medium Deep Survey Team at JHU, especially Eric J. Ostrander 
for his efforts on retrieving and reducing the archival HST data. We are grateful to Emilio Falco 
for providing the lens magnitudes of 0142-100 system. We also thank Chris Kochanek, Howard 
K. C. Yee, Joel Primack, and Stefano Casertano for useful discussions and communications, and 
Mark Subbarao and the anonymous referee for helpful comments on the manuscript. 



TABLE 1 
List of Gravitational Lenses 



Name 


rriL (total) 


zl 


zs 


O.r.ti") 


References 


HST14176 + 5226 


1=19.71 ± 0.05 


0.81 


3.4(') 


1.51 


Ratnatunga et al. (1995), Crampton et al. (1996) 


HST12531-2914 


1=21.82 ± 0.05 


0.7 ± 0.1 


< 5 


0.65 


Ratnatunga et al. (1995) 


PG1115 + 080 


R=18.36 ± 0.3 


0.29 


1.72 


1.10 


Weymann et al. (1980), Kristian et al. (1993) 


MG1654+1346 


R=18.4 ± 0.3 


0.25 


1.74 


1.05 


Langston et al. (1989) 


CLASS1608+656 


R=19.7 ± 0.3 (K~ 16) 


0.63 


1.39 


1.05 


Myer et al. (1995a), Fassnacht et al. (1996) 


0142-100 


R=19.36 ± 0.10 


0.49 


2.72 


1.10 


Surdej et al. (1987), Falco (1995) 


MG0414+0534 


1=21.22 ± 0.15 


1.2 ± 0.4 


2.64 


1.05 


Hewett et al. (1992), Lawrence et al. (1995), Schechter & Moore (1992) 


B0218+357* 


R=20.0 ± 0.3 


0.68 


0.94 


0.35 


Patnaik et al. (1993) 



* Spiral galaxy lens 



- 11 - 



REFERENCES 

Alock, C, & Anderson, N. 1986, ApJ, 302, 43 

Barrientos, L. F., Schade, D., & Lopez-Cruz, O. 1996, ApJ, 460, L89 
Baugh, C. M., Cole, S., & Frenk, C. S. 1996, MNRAS, submitted 
Bender, R., Ziegler, B., & Bruzual, G., ApJ, 463, L51 
Breimer, T. G., & Sanders, R. H. 1993, A&A, 274, 96 
Bruzual, G. A., & Chariot, S. 1993, ApJ, 405, 538 

Burke, B. F. 1990, in Lecture Notes in Physics, 360, Gravitational Lensing, ed. Y. Meillier, B. 
Fort, h G. Soucail, (Berlin: Springer), 127 

Crampton, D., Le Fevre, O., Hammer, F., h Lilly, S. J. 1996, A&A, in press 

Fabcr, S., k Jackson, R. 1976, ApJ, 204, 668 

Falco, E. 1995, private communication 

Fassnacht, C. D., et al. 1996, ApJ, in press 

Franx, M. 1993, in Galactic Bulges, ed. H. Dejonghe k, H. J. Habing (Dordrecht: Kluwer) 243 

Freedman W. L. et al. 1994, Nature, 371, 757 

Fukugita, M. & Peebles, P. J. E. 1995, Princeton preprint 

Fukugita, M., Futamase, T., Kasai, M., & Turner, E. L. 1992, ApJ, 393, 3 

Fukugita, M., Futamase, T., & Kasai, M. 1990, MNRAS, 246, 24P 

Gott, J. E. IIL, Park, M-G., & Lee, H. M. 1989, ApJ, 338, 1 

Helbig, P., & Kayser, R. 1996, A&A, in press 

Hammer et al. 1995, A&A, 298, 737 

Hewitt et al. 1992, AJ, 104, 968 

Im, M., Griffiths, R. E., Ratnatunga, K. U., & Sarajedini, V. L. 1996, ApJ, 461, L79 
Im, M., Griffiths, R. E., & Ratnatunga, K. U. 1997, ApJ, submitted 
Impey, C. D. et al. 1996, ApJ, in press 

Kauffmann, G., Chariot, S., & White, S. D. M. 1996, MNRAS, submitted 

Keeton II, C. R., k Kochanek 1995, in lAU Symposium 173: Astrophysical Applications of 
Gravitational Lensing, ed. C.S. Kochanek k J.N. Hewett, (Dordrecht:Kluwer) 419 

Kochanek, C. S. 1995, ApJ, submitted 

Kochanek, C. S. 1994, ApJ, 436, 56 

Kochanek, C. S. 1993, ApJ, 419, 12 



- 12 - 



Kochanek, C. S. 1992, ApJ, 384, 1 

Krauss, L. M., & White, M. 1992, ApJ, 394, 385 

Kristian, J. et al. 1993, AJ, 106, 1330 

Langston et al. 1989, AJ, 97, 1283 

Lauer, T. 1988, ApJ, 325, 49 

Lawrence, C. R., Cohen, J. G., & Oke, J. B. 1995, AJ, 110, 2583 

Loveday, J., Peterson, B. A., Efstathious, G., & Maddox, S. J. 1992, ApJ, 390, 338 

Maoz, D., k Rix, H.-W. 1993, ApJ, 416, 425 

Marzke et al. 1994, AJ, 108, 437 

Miralda-Escude, J. 1991, ApJ, 370, 1 

Myers, S. T. et al. 1995, ApJ, 447, L5 

Ostriker, J. P., & Steinhardt, P. J. 1995, preprint 

Paczynski, B., & Gorski 1981, ApJ, 248, LlOl 

Pahre, M. A., Djorgovski, S. G., & de Carvalho, R. R. 1996, ApJ, 456, L79 
Patnaik, A. R., et al. 1993, MNRAS, 261, 435 

Ratnatunga, K. U., Ostrander, E. J., Griffiths, R. E., & Im, M. 1995, ApJ, 453, L5 
Ratnatunga, K. U., et al. 1996, in preparation 

Rees, M. 1996, in the Cosmological Constant and the Evolution of the Universe, Proc. First 
RESCUE Symposium, ed. M. Fukugita (Universal Academy Press, Tokyo) 1 

Schechter, P. L., & Moore, C. B. 1992, AJ, 105, 1 

Surdej, J. & Soucail, G. 1994, in Gravitational Lenses in the Universe, ed. J. Surdej, D. 
Praipont-Caro, E. Gosset, S. Refsdal, & M. Remy (Liege, Universite de Liege) 153 

Surdej, J. et al. 1987, Nature, 329, 695 

Turner, E. L. 1990, ApJ, 365, L43 

Turner, E. L., Ostriker, J. P., & Gott, J. R. IIL 1984, ApJ, 284, 1 

Wallington, S., & Narayan, R. 1993, ApJ, 403, 517 

Weymann, R. J. et al. 1980, Nature, 285, 641 

Yee, H. K. C. 1995, private communication 

Yee, H. K. C, & Ellingson 1994, AJ, 107, 28 

de Zeeuw, T., & Pranx, M. 1991, ARA&A, 29, 239 



This preprint was prepared with the AAS lAT^^X macros v3.0. 



